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Abstract 

The paper concerns the time-dependent behavior of electroactive polymers (EAP) and their 
use in advanced intelligent structures for space exploration. Innovative actuator design for low 
weight and low power valves required in small plants planned for use on the moon for chemical 
analysis is discussed. It is shown that in-depth understanding of cyclic loading effects observed 
through accelerated creep rates due to creep-fatigue interaction in polymers is critical in terms of 
proper functioning of EAP based actuator devices. In the paper, an overview of experimental 
results concerning the creep properties and cyclic creep response of a thin film piezoelectric 
polymer polyvinylidene fluoride (PVDF) is presented. The development of a constitutive 
creep-fatigue interaction model to predict the durability and service life of electroactive 
polymers is discussed. A novel method is proposed to predict damage accumulation and 
fatigue life of polymers under cyclic loading conditions in the presence of creep. The study 
provides a basis for ongoing research initiatives at the NASA Kennedy Space Center in the 
pursuit of new technologies using EAP as active elements for lunar exploration systems. 
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Introduction 


In the past decades electroactive polymers (EAP) have been increasingly used in many 
technological fields including aerospace, transportation, telecommunication, photonics, and 
bioengineering due to their unique ability to sense and react to external stimuli in the form of 
electrical signals or mechanical deformations. Polymer based transducers, sensors and actuators 
have shaped an impressive array of advanced technologies including acoustic microscopy, active 
vibration control, damage detection in fibrous composites, optoelectronics, and marine fouling 
prevention. Increasingly, electroactive polymer systems have been used for space exploration in 
weight sensitive applications such as shape and position control of compliant structures, smart 
skins, solar sails, deployable membrane mirrors, atmospheric balloons, antennae, reflectors, 
flexible robotic manipulators for grasping and locomotion, pumps, valves, relays and other 
applications which have traditionally relied on electromagnetic components. The use of 
electroactive polymers provides numerous advantages in terms of reduced weight, large motions 
without gearing, and built in position sensing capabilities [1,2]. 

The National Aeronautics and Space Administration is currently pursuing a lunar 
exploration program that will benefit from innovations in system components designed to reduce 
weight and power consumption without sacrificing reliability. For example, low weight and low 
power fluid handling valves will be required in small chemical plants planned for use on the 
moon for chemical analysis. Available valves based on solenoids can produce millions of 
opening and closing cycles so that valve seat design rather than the electromagnets becomes the 
limiting life design factor. In these applications, electroactive polymer actuators are considered 
to be a logical substitute for current electromagnetic devices. 

The effective use of electroactive polymers as active elements in actuator design requires 
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precise understanding of their properties. Many important issues in this regard arise in 
connection with material nonlinearities, time-dependent effects, and energy losses that 
characterize the electromechanical response of EAP [3-7]. It is well known that under sustained 
loading conditions polymers exhibit creep at room temperature with accelerated creep rates at 
elevated temperatures. The creep behavior of polymers in response to static loads changes 
dramatically in response to cyclic load histories involving stress reversals or cyclic temperatures. 
Under such conditions polymers develop much higher creep rates than those observed in static 
experiments. It is considered that cyclic creep acceleration in polymers is caused by material 
damage due to fatigue which results in creep-fatigue interaction [8, 9]. 

Enhanced understanding and accurate prediction of creep-fatigue interaction effects on 
the performance and longevity of EAP devices is critical. Since such devices are designed to 
operate in dynamic environments, accelerated material degradation processes are likely to impair 
the functional performance of EAP systems, similarly to the effects of microcracking that tend to 
decrease the accuracy of piezoelectric sensors and actuators by weakening the electromechanical 
couplings in material properties [10]. 

Experimental Results 

A representative group of EAP comprises piezoelectric polymer systems including 
polyvinylidene fluoride (PVDF), PVDF co-polymers and piezoelectric polymer based 
composites. At present, this group has reached a greater degree of maturity in terms of 
technological applications as compared with other electroactive polymers. Essentially, PVDF 
and PVDF co-polymers represent the principal commercially available high-performance 
piezoelectric polymers currently in use. These materials are characterized by many attractive 
properties such as stable response characteristics in a wide frequency range of up to 10 9 Hz, low 
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acoustic impedance, high degree of resistance to impact, resistance to moisture absorption, and 
marginal sensitivity to ultraviolet and nuclear radiation [4-6]. 

Strong piezoelectric effects in PVDF are obtained through a technological process that 
involves stretching and polling of extruded thin sheets of the polymer. Stretching provides an 
alignment of molecular chains in the stretch direction. An applied electric field of 100 kV/mm at 
an elevated temperature of about 103°C causes permanent polarization which is maintained after 
the material is cooled to room temperature. These processing and fabrication conditions play a 
decisive role in defining the electromechanical properties of piezoelectric polymers. 

The time-dependent mechanical properties of PVDF have been studied based on a series of 
quasi-static creep tests and dynamic mechanical tests of 28 pm thick commercially produced 
PVDF thin films. In the experiments, PVDF samples have been tested in two in-plane material 
directions, i.e., parallel (direction 1) and perpendicular (direction 2) to the aligned molecular chains 
of the polymer. Creep experiments have been performed at 10 different stress levels under 
sustained loading conditions at room temperature. Strain measurements have been taken by a linear 
variable differential transformer (LVDT). The creep behavior of PVDF has been also studied 
using the method of dynamic mechanical testing and analysis. The results of these experiments are 
illustrated in Figs. 1 and 2 that demonstrate the creep response of PVDF in terms of strain e (%) as 
a function of time t [1 1-12]. 

It has been determined that at stresses a below 57% of the yield stress a yi = 45 MPa in 
the material direction 1, and below 76% of the yield stress a y 2 = 39 MPa in the direction 2 the 
creep response of PVDF is linear such that the strain histories e(t) obtained at various constant 
stresses and normalized based on the respective stress can be represented by a single curve 
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known as creep compliance. The creep compliances for PVDF, C„(t) = e(t)/a (n= 1, 2), for 
material directions 1 and 2 are shown in Figs. 3 and 4, respectively. 



Fig. 1. Creep Strain of PVDF (direction 1). 



Fig. 2. Creep Strain of PVDF (direction 2). 
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Fig. 3. Creep Compliance of PVDF (direction 1). 
C 2 (t) (1/Pa) 



Fig. 4. Creep Compliance of PVDF (direction 2). 

Cyclic creep response of PVDF has been studied in the stretch direction of the polymer 
(direction 1) because of higher creep rates as compared with the direction 2. The experimental 
investigation involved a series of tests under the conditions of superimposed static and cyclic 
loads, i.e., a(t) = a m + <r a sin cot applied in asymmetric tension-tension mode. The applied cyclic 
amplitude, a a was much smaller than the mean stress <7 m , i.e. a a « a m . The results were obtained 
in terms of creep strains e(t) and interpreted in terms of normalized creep strain e(t)/a m . At least three - 


6 





five identical experiments were repeated under the same loading conditions to ensure reproducible 
results [9, 13]. 

This study has demonstrated that under cyclic loading conditions creep rates of PVDF 
accelerated with an increase of the amplitude and frequency of the cyclic load. These effects are 
illustrated in Fig. 5 which provides a comparison between cyclic creep strains e(t) normalized by 
the mean stress a m and the respective creep compliances obtained under static loading conditions. 
Cyclic creep of PVDF depicted by curve 3 has been obtained for the following loading conditions: 
o m = (0.45)oy, c a = (0.01 4 )oy, and co = 10 Hz. The curve 4 was obtained for o m = (0.34)oy, 
a a - (0.068) ay, and to = 5 Hz. The curves 1 and 2 represent static creep compliances obtained at 
a = (0. 34 )oy and a = (0. 45)oy, respectively. 

(e/o m ) xlO lo (l/Pa) 



Fig. 5. Cyclic and Static Creep Responses of PVDF. 


These and other results reported in [9, 13] demonstrate that piezoelectric polymers 
undergo measurable creep acceleration when subjected to superimposed static and cyclic loads. 
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Consistently, cyclic creep rates tended to increase with an increase of frequencies or amplitudes. 
This phenomenon can be explained by the presence of several interactive mechanisms: creep, 
hysteretic heating, and fatigue damage. The contribution of each individual mechanism depends on 
the loading and temperature conditions. In particular, experiments indicate that fatigue failure of 
polymers is dominated by hysteretic heating at higher stress levels and frequencies, whereas 
damage accumulation processes become critical at lower stresses and frequencies [8, 14], In 
general, cyclic damage evolution in polymers involves several consecutive stages, the formation 
of crazes, craze growth, crack nucleation and crack propagation [14]. Even before the appearance 
of crazes, polymers subjected to cyclic loading regimes tend to undergo measurable changes of 
their density, shear modulus, and internal damping. The study reported in [15] has shown direct 
correlation between the sensitivity of polymers to cyclic frequency effects and damage evolution 
in terms of nucleation and propagation of crazes. These observations, however, have not been 
consistently characterized or quantified. To date, the creep-fatigue interaction effects in the 
behavior of electroactive polymers remain unexplored. 

Modeling of Static Creep 

As discussed in the foregoing section, within the range of stresses not exceeding 25 MPa the 
creep properties of PVDF can be represented by creep compliances shown in Figs. 3 and 4. 
Respectively, the polymer can be characterized by the constitutive equations of linear 
viscoelasticity based on the Boltzmann superposition principle which states that since strain is a 
linear function of stress the total effect of applying several stresses is the sum of the effect of 
applying each one separately. The application of this principle leads to the stress-strain relation 
in the form [16, 17] 
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e(t) = a(0)C n (t) + Jc n (t-T)^dx, 


dx 


0 ) 

where the stress a and strain e are functions of time t, and C n (n = 1, 2) denote the creep 
compliances in the respective material direction of PVDF. 

The creep compliances C n (t) of PVDF can be represented analytically in the form 

Cn (0 = + b n (tn ), n— l, 2 


( 2 ) 


where the coefficients a n , b n and a n are defined as follows: 

Direction 1: a, = 3.206xl0' 10 ; b, = 5.018xl0' u ; ai = 0.107 (3) 

Direction 2: a 2 = 3.514xlO' 10 ; b 2 = 0.11 1x10'"; <x 2 = 0.085 (4) 

Substitution of Eq. (2) into Eq. (1) leads to an integral equation which can be solved 
using the Laplace transformation method in terms of the stress function a(t) in the form 


*«) 1} S (Wa„)r(a n )(t-x) a “ rV . 

<*(<)=- _\r~/ s <^ T ;n=l,2 

n n n 1 


(t-x)r(na n ) 

where T (a n ) and T (na n ) are gamma functions, and the coefficients a n , b n , and On (n = 1,2) are 
defined by Eqs. (3 and 4). 

In the case of constant strain e(t) = e 0 = const., Eq. (5) is of the form 


(5) 




1 \ 2 (b n a n /a„ )r(a_ )(t - x)“ n 1 

o(t)/e 0 = — [1- J ° ~ <1t]= — [1 - R„(t)] ; n =1> 2 


o n = 1 


(t-x)T(na n ) 


( 6 ) 


where 


! ? On«„/a„)r(a„)(t-x) a " 

Rn(t)=J . r . dt ;n = l,2 

q n- 1 (t-x)r(na n ) 


( 7 ) 
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denote the relaxation functions of PVDF in both in-plane material directions. The set of 
Eqs. (1-7) completes the characterization of PVDF in the linear viscoelastic range of creep 
properties under static loading conditions. 

Modeling of Creep-Fatigue Interaction 

To characterize creep-fatigue interaction effects in PVDF, a constitutive model has been 
developed based on the principles of linear viscoelasticity and continuum damage mechanics 
[16-19]. The concept of a fictitious undamaged continuum with effective material characteristics 
has been used to represent fatigue damage in the material. Specifically, damage evolution in 
PVDF has been characterized by a damage function D which depends on the number of cycles 

N, i.e., D = D (N). The damage function is defined within the limits 0 < D < D*, where the low 

* 

limit represents undamaged material, D(0) = 0, and D characterizes the degree of damage at 

which the response of the material cannot be treated anymore as linearly viscoelastic. The value 
* 

ofD must be determined experimentally. 

Based on these assumptions, the properties of a fatigue damaged linear viscoelastic 
material can be defined by the following constitutive equation 

e D (t, N) = Ceff (t, N)g( 0) + } C eff (t - x, N)-j-dx (8) 

t dx 

where 

C e fi(t, N) = C(t)/[ 1 - D(N)] (9) 

denotes the effective creep compliance of the material with damage and C(t) is the creep 
compliance of undamaged material. Note that, at N = 0, C e ff(t, 0) = C(t). In Eq. (8), e D (t, N)) 
denotes the strain in the damaged material produced in response to the stress history a(t). 
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It follows from Eq. (9) that the damage function D(N) can be represented in the form 


D(N)= 1 - 


C(t) 


( 10 ) 


Based on this representation an experimental protocol can be developed to 
characterize the damage function D(N) since the creep compliances C e ff(t, N) and C(t) can be 
determined experimentally. Specifically, the creep compliance C(t) can be obtained from 
static creep experiments at a constant stress ct 0 by normalizing creep deformations based on 
the applied stress. The same method can be applied to a damaged material after N number of 
cycles, providing the effective creep compliance C e ff (t, N). Note that this approach is based on 
the assumption that the damaged material is linearly viscoelastic. 

Since the creep compliances C(t) and C e ff(t, N) are obtained by normalizing the 
experimentally determined strains in the damaged and undamaged materials, 8o(t, N) and e(t), 
respectively, such that 

Ceff (t, N) = e D (t, N)/a 0 and C (t) = e(t)/a 0 (11) 

substitution of Eqs. (1 1) into Eq. (10) provides an alternative equation for computing the 
damage function 

D(N)=l--7%r (12) 

The proposed creep-fatigue interaction model has been validated experimentally. PVDF 
samples were prepared from commercially available thin sheets of the material with deposited 
silver electrode layers on both surfaces of the polymer. The total thickness of the samples 
including the thickness of silver layers was 46 pm, with the thickness of the PVDF layer of 28 
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pm. The samples had in-plane dimensions 22 mm x 1 .5 mm, and were tested in the direction of 
the aligned molecular chains of the polymer (direction 1). Tests were performed using a TA 
Instruments dynamic mechanical analyzer DMA 2980. 

As the first step, the creep compliance C(t) = e(t)/a 0 of undamaged material was . 
determined from creep tests under sustained loading conditions at the stress level 
(j 0 = 10 MPa, which is below the viscoelastic linearity limit. The results of these tests are 
shown by the curve 1 in Fig. 6. 

Creep Compliance (1/Pa)x10' 6 



Fig. 6. Creep ofPVDF (1: Static; 2: Cyclic ; 3: Two-Stage Stress History). 

The second series of experiments was conducted under the conditions of superimposed 
static stress cj 0 = 10 MPa and cyclic stress with the amplitude a a = 4.5 MPa and frequency 
<d = 20 Hz. Note that the total maximum stress Omax = 14.5 MPa did not exceed the viscoelastic 
linearity limit ofPVDF determined as 25.65 MPa. The respective cyclic strains were measured 
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and normalized with respect to a 0 as represented by the curve 2 in Fig. 6. As expected, at this 
stage, the material demonstrated visible creep acceleration as compared with static creep. 

The third series of experiments consisted of two types of tests. Firstly, cyclic damage was 
generated for over six hours under the same stress conditions as applied at the second stage of 
the program. After that period only the static load was maintained. The respective effective creep 
compliances C e ff (t) was determined by normalizing the measured creep strain of the damaged 
material by the respective static stress. The results of these experiments are illustrated by the 
curve 3 in Fig. 6. 

It is clear that under static loading conditions the damage function Di = 0. Based on the 
experimental data obtained from the second and third series of tests, the respective damage 
functions D2 and D3 shown in Fig. 7 were determined using Eq. (10). 



Cycles 

Fig. 7. Damage Characterization of PVDF. 

It follows from these diagrams that cyclic damage in PVDF initiates immediately upon 
the load application. Both functions D2 and D3 increase with the number of cycles N which 
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characterizes the rate of damage evolution. The damage function D 3 remains constant, which 
indicates no further damage development in the material in the absence of the cyclic load. 

Challenges 

Integration of electroactive polymers into structural design provides the capability of 
controlling the mechanical characteristics of structures in terms of stiffness or damping, or 
modify the structural response in terms of position or velocity. This type of built-in structural 
intelligence has been particularly effective in space applications. In particular, the use of EAP in 
machine components and devices for lunar exploration has been actively pursued at the NASA 
Kennedy Space Center. In November of 2008, NASA performed a series of experiments at the 
lunar analog site on Mauna Kea, Hawaii, chiefly to demonstrate the technologies capable of 
analyzing and extracting oxygen from lunar soil. Experiments like the Regolith Volatiles 
Characterization (RVC) and the Lunar Water Resource Demonstrator (LWRD) were essentially 
miniaturized chemical plants for detecting soil volatiles and analyzing water [20]. The plants 
contained rotary valves, solenoid valves and pumps for routing gases. Ambitious goals have been 
set for the design of these devices in terms of mass, thermal, and power requirements. In this 
regard, the existing technologies appear ineffective. It is expected that a promising alternative to 
traditional design will be provided by electroactive polymer systems that will enable new 
technological developments for lunar exploration. 

It is clear that, in practice, successful implementation of EAP systems directly depends 
on the degree of understanding of their behavior and properties. To date, considerable progress 
has been made in this subject area [2, 4]. However, many aspects of the performance and 
properties of electroactive polymers remain unexplored. In particular, challenges arise due to 
the sensitivity of polymers to fabrication and temperature conditions, time-dependent effects, 
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and material nonlinearities. An immediately obvious challenge to fielding electroactive 
polymers is a lack of understanding of their failure modes, reliability, and the change of 
electromechanical properties depending on time and temperature. Questions as to the nature of 
long term creep and fatigue-fatigue interaction effects in polymers must be understood and 
quantified before EAP can be widely accepted as viable competitors to the electromagnetic 
actuator technologies currently in use. 


Work-in-Progress 

To investigate the potential use of EAP in lunar exploration, the Kennedy Space Center is 
sponsoring research on the creation of simulation models to assist in the design of EAP based 
active control systems. To illustrate the importance of understanding long term creep and fatigue 
effects in an EAP element, consider the design of a valve currently under investigation at 
Kennedy Space Center. A simplified valve schematic and the corresponding valve seal free body 
diagram demonstrating the acting forces are shown in Fig. 8. 


W = membrane width (constant) 

◄ ► 




F v viscoelastic forces 
F f frictional forces 
F r valve seat reactive force 
F g force of gas pressure 

F s spring force 

F m Maxwell (electrostatic) force 


Fig. 8. Valve Schematic and Free Body Diagram. 
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In this normally open design, a valve seal moves downward on application of voltage to 
mate with a seat and cut off the flow of gas. When no voltage is applied, the seal is at xi=0, and 
the spring forces F s and membrane viscoelastic forces F v are in equilibrium. As Maxwell force 
F m increases with increasing voltage, the seal moves toward the seat to a new equilibrium 
position. At sufficient voltage the Maxwell force overcomes viscoelastic, frictional, and gas 
pressure forces so that the seal contacts the seat resulting in a reactive force at the seal-seat 
interface 

F r = F s + F m - F v - F g - F f (13) 

The reactive force F r , equivalent to the blocking force in [21], for a given applied voltage 
can be increased by designing the seat as close as possible to the seal. However, due to the long 
term effects of creep and fatigue the initial zero-voltage equilibrium position of the seal will drift 
towards the seat so that designing with the initial seal position too close to the seat would cause 
the valve to fail prematurely. It is clearly important to be able to predict these time-dependent 
effects under the given loading and temperature conditions in order to accurately design the seal 
to seat gap to avoid gas leakage and increase the reliability of the valve. 

Potentially, an electroactive EAP membrane can be used as a sensor and modeled as a 
capacitor with the capacitance being a function of the actuator geometry and the dielectric 
properties of the polymer [22]. As the membrane extends, the capacitance will increase in a 
predictable and repeatable manner. With proper design of circuitry, an EAP valve can exploit 
this built-in displacement sensor without the need for a separate sensing device. It should be 
noted that since highly piezoelectric materials like PVDF have the additional property of 
generating voltages when deformed, this property could also be used to measure actuator 
position. However for static or very slow frequency measurements, determination of position 
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from measurement of capacitance may prove to be more accurate and will be considered here for 
the fault detection algorithm. 

Long-term drift can be measured and compared to the predicted creep models. 
Specifically, and estimated damage function, D(n ) , can be updated using the empirically 
determined model and a count of total valve cycles, n. When the valve is closed, the constant 
strain is known, and when the valve is open, the stress is known, based on the known force Fs 
and the respective position of the valve. Using the material characterization provided by the 
constitutive Eqs. (7) and (8) together with the measured opening and closing times and positions, 
it will be possible to simulate an estimated strain, s(n) estimate the accumulated damage and 
compare this estimate to the measured strain giving the rule 

|e(n) - e(n)| < k (14) 

where k denotes a certain limit based on the accuracy of the estimates, and violation of the rule 
signifies a valve fault. In practice, accurate simulation of e(n) may be computationally 
intensive. Instead, an approximation consisting of a series of exponential terms will be sought 
based on the analysis of equations (7) and (8). It is hypothesized that this estimator could be 
designed with a low order digital filter and would not be computationally intensive. 

Conclusions 

Experimental results show that the mechanical properties of electroactive polymers are 
time-dependent. Within certain limits creep deformations of EAP measured at room temperature 
can be accurately predicted using the constitutive equations of linear viscoelasticity. The 
dynamic response of EAP subjected to superimposed static and cyclic loads is characterized by 
accelerated creep rates due to the effects of creep-fatigue interaction. Cyclic creep acceleration 
has been observed even in the range of stresses well below the viscoelastic linearity limit. It is 
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clear that the cyclic response of the EAP is essentially nonlinear, since it does not represent a 
simple superposition of the responses to static and fully reversed cyclic loads applied separately. 

To characterize cyclic damage evolution in PVDF, a constitutive material model has been 
developed based on the principles of linear viscoelasticity and continuum damage mechanics. 

The model has been validated experimentally. Fatigue induced material degradation has been 
characterized by experimentally determined damage function. This approach provides a reliable 
predictive capability for assessing the long term integrity and functionality of EAP. 

The study presented in this paper provides a basis for ongoing research initiatives at the 
NASA Kennedy Space Center in the pursuit of new technological developments using EAP as 
active elements for lunar exploration systems. 

References 

[1] Ihlefeld, C. M. and Qu, Z., “A Dialectric Electroactive Polymer Generator-Actuator Model: 
Modeling, Identification, and Dynamic Simulation,” Proceedings of SPIE "Electroactive 
Polymer Actuators and Devices (EAPAD) 2008, ” (Y. Bar-Cohen, editor), San Diego, CA, March 
10-13, SPIE, 2008, Paper 6927-OR. 

[2] Vinogradov, A. M., “Accomplishments and Future Trends in the Field of Electroactive 
Polymers” (Invited), Proceedings of SPIE "Electroactive Polymer Actuators and Devices (EAPAD) 
2008, ” (Y. Bar-Cohen, editor), San Diego, CA, March 10-13, 2008, SPIE, Paper 6927- OM. 

[3] Srinivasan, A.V., and McFarland, D. M., 2001, “Smart Structures. Analysis and Design,” 
Cambridge University Press, Cambridge, UK. 

[4] Vinogradov, A. M., “Piezoelectricity in Polymers,” Encyclopedia of Smart Materials, (M. 
Schwartz, editor) Vol. 2, J. Wiley 8c Sons, 2002, pp. 780-792. 


18 



[5] Zhang, Q.M., V. Bharti, V., and Kavamos, G., “Poly(vinylidene fluoride) (PVDF) and its 
Copolymers,” Encyclopedia of Smart Materials, (M. Schwartz, editor) Vol. 2, J. Wiley & Sons, 
NY, 2002, pp. 807-825. 

[6] Harrison, J.S. and Ouriaies, Z. “Polymers, Piezoelectric,” Encyclopedia of Smart Materials, 
(M. Schwartz, editor) Vol. 2, J. Wiley & Sons, NY, 2002, pp. 860-873. 

[7] Vinogradov, A. M., Schmidt, V. H., Tuthill, G. F., and Bohannan, G. W., “Damping and 
Electromechanical Energy Losses in the Piezoelectric Polymer PVDF,” Mechanics of Materials, 
Vol. 36, 2004, pp. 1007- 1016. 

[8] Lesser, A. J., “Changes in Mechanical Behavior During Fatigue of Semicrystalline 
Thermoplastics,” Journal of Applied Polymer Science, Vol. 58, No. 5, 1995, pp. 869 -879. 

[9] Vinogradov, A. M. and Schumacher, S.C., “Cyclic Creep of Polymers and Polymer Matrix 
Composites,” Mechanics of Composite Materials, No. 1, 2001, pp. 29-34. 

[10] Sevostianov, I., and Kachanov, M., “Microcracking in Piezoelectrics Weakens the Electro- 
Mechanical Coupling and Changes its Directionality,” International Journal of Fracture, Vol. 
101, 2000, pp.L-3-L-8. 

[11] Holloway, F., “Material Characterization of Poly(vinylidene Fluoride): A Thin Film 
Piezoelectric Polymer,” M.S. Dissertation, Montana State University-Bozeman, 1997. 

[12] Vinogradov, A. M., and Holloway, F., “Electro-Mechanical Properties of the 
Piezoelectric Polymer PVDF,” Ferroelectrics, No. 226, 1999, pp. 169-171. 

[13] Schumacher S.C. “Parametric Study of Cyclic Loading Effects on the Creep Behavior of 
Polymers and Polymer Based Composites,” MS Dissertation, MSU, 2000. 

[14] Hertzberg, R. W., and Manson J. A., ‘Tatigue of Engineering Plastics,” Academic Press, 1980. 


19 



[15] Hertzberg R.W., J. A. Manson, J. A., and M. Skibo, “Frequency Sensitivity of Fatigue 
Processes in Polymeric Solids,” Polymer Eng. Sci., Vol. 15,No. 4, 1975, pp. 252-260. 

[16] Ferry, J. D., “Viscoelastic Properties of Polymers,” 3-rd ed., J. Wiley & Sons, 1980. 

[17] Wineman, A. S. and Rajagopal, K. R., “Mechanical Response of Polymers,” Cambridge 
Univ. Press, 2000. 

[18] Krajcinovic, D., “Damage Mechanics,” Elsevier, 1996. 

[19] Skrzypek, J., and A. Ganczarski, A., “Modeling of Material Damage and Failure of 
Structures,” Springer, 1 999. 

[20] Muscatello, A., Captain, J., Lueck, D., Quinn, J., Gibson, T., and Weis, K., “Lunar Water 
Resource Demonstrator,” Proceedings of 47 th AIAA Aerospace Sciences Meeting, Orlando, FL, 
January 5-8, 2009, Paper AIAA-2009-1202, 2009. 

[21] Kofod, G., Sommer-Larsen, P., Kombluh, R., Pelrine, R., “Actuation Response of 
Polyacrylate Dielectric Elastomers,” Journal of Intelligent Material Systems and Structures, 

Vol. 14, Dec. 2003, pp.787-793 

[22] Pelrine, R., Kombluh, R., Eckerle, J., Jeuck, P., Oh, S., Pei, Q. Stanford, S., “Dielectric 
Elastomers: Generator Mode Fundamentals and Applications,” Proceedings of the SPIE “ Smart 
Structures and Materials 2001: Electroactive Polymer Actuators and Devices, ” (Y. Bar-Cohen, 
editor) Vol. 4329, SPIE, 2001, pp. 148-156. 


20 



